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Abstract 

The dynamics of linear stochastic growth equations on growing substrates is studied. The sub- 
strate is assumed to grow in time following the power law t 7 , where the growth index 7 is an 
arbitrary positive number. Two different regimes are clearly identified: for small 7 the interface 
becomes correlated, and the dynamics is dominated by diffusion; for large 7 the interface stays 
uncorrelated, and the dynamics is dominated by dilution. In this second regime, for short time 
intervals and spatial scales the critical exponents corresponding to the non-growing substrate sit- 
uation are recovered. For long time differences or large spatial scales the situation is different. 
Large spatial scales show the uncorrelated character of the growing interface. Long time intervals 
are studied by means of the auto-correlation and persistence exponents. It becomes apparent that 
dilution is the mechanism by which correlations are propagated in this second case. 

PACS numbers: 68.35.Ct,05.40.-a,64.60.Ht 
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I. INTRODUCTION 



Fluctuating interfaces have been the object of study of many different works over the 
last decades. Together with the possible technological applications that their understanding 
may bring, as for instance in thin film industry, there is a genuine theoretical interest in 
unveiling their dynamical properties. This is so to the extend that the Kardar-Parisi- Zhang 
(KPZ) equation [if], one of the most influential models for surface growth, is being currently 
considered as a prototypical model of nonequilibrium dynamics. 

Usually, stochastic equations modeling surface growth have been studied in static do- 
mains. On the other hand, some types of growing interfaces, as for instance radial ones, 
present a domain size that grows over time. The direct study of radial interfaces is com- 
plicated by nonlinear effects, including the possibility of instabilities affecting the radial 
symmetry j^, 3, 4|. This suggests studying first the dynamics of linear stochastic growth 
equations on growing domains. Once the effect of substrate growth on the interface dynamics 
is understood, one could move to the more complicated case of radial growth. 

Of course, considering linear growth equations has a limited applicability to real physical 
systems, as important nonlinearities are being neglected. Still, the detailed analysis of linear 
stochastic growth equations has revealed important physical properties of rough surfaces . 
Additionally, the study of linear equations has served as a basis for the approach to the more 
complicated nonlinear ones [6j]. In this sense, we expect that the results presented here might 
bring useful insights into the dynamics of nonlinear equations on growing domains. 

Stochastic growth equations for radial interfaces have been previously considered in the 
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literature [2|, |3|, U, [7|, [8||. In [8], radial interfaces with a linearly in time growing domain 
are studied, and several dynamical quantities are calculated and compared to the classical 
values. Under the approximations made in this work, the only genuine radial effect that 
is being considered is domain growth, while the possible nonlinear effects are disregarded. 
This allows a perfect comparison among the results presented there and the ones that we 
will introduce here. We will also examine how decorrelation might appear in the growing 
interface and what is the resulting large scale structure 4J, how the classical values of the 
critical exponents are recovered and in what precise limits this occurs 1(3 ] . 

The goal of this work is to put previous studies considering linearly in time growing 
domains in a broader context. This will be achieved by proposing an arbitrary power law 
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growth model for the domain. We will focus on rough interfaces, i. e., models for which the 
growth exponent (3 is strictly positive. The paper is organized as follows: in section [Til we 
describe the phenomenology of domain growth. In sec. II I II we focus on the nonequilibrium 
dynamics of the Edwards- Wilkinson (EW) equation, and in sec. HVI we extend these results to 
the general linear Langevin dynamics. In sec. |V]the temporal correlations of the fluctuating 
interface are studied, and in sec. |VT] its persistence properties. In sec. IVHI the connection 
to radial growth is investigated, and the conclusions of this work are drawn in sec. IVIIIl 



II. GROWING DOMAINS 



In order to study the dynamics of stochastic growth equations on growing domains we 
begin considering the EW equation which reads 

d t h = DV 2 h + F + ay,t), (1) 

where £(y,i) is a zero-mean Gaussian white noise which correlation is 

(e(y,W,0) = ^(y-y / )^-t / ), (2) 

D is the diffusion constant, F the constant deposition rate and e the noise intensity, all these 
parameters being positive. To derive the EW equation on a growing domain we will follow 
the theory introduced in 121 ]. focused on reaction diffusion dynamics on uniformly growing 
domains. We start considering the conservation law in integral form 

/ h(y,t)dy= [ [_V-j + ^(y,t)]rfy, (3) 
at J s t J s t 

where St is the uniformly growing domain, j = -DV/j is the current generated by diffusion, 
and jF(y,t) = F + £(y, t) is the EW growth mechanism. By applying Reynolds transport 
theorem we find 



d 
dt 



[ h{y,t)dy= [ [d t h + V-(vh)]dy, (4) 
Js t Js t 

where v(y, t) denotes the flow velocity generated by the growing domain. Valid as it is for 
any domain, the integral conservation law may be expressed in the local form 

d t h + V- (vh) = DV 2 h + F(y,t). (5) 

In this equation we readily identify two new terms, the advection one v- V/t, and the dilution 
one hV ■ v. For every y G St, that has evolved from y G S to , we find v(y,t) = dy/dt. Let 



us now concentrate in one- dimensional substrates and then move to higher dimensionalities. 
In this case uniform growth translates into y = g(t)yo, where g(t) is a temporal function 
such that g(t ) = 1. This yields v = yg/g, and thus 

dth + ft (yd y h + h) = Dd 2 y h + F + £(</, t). (6) 

For a one-dimensional substrate (0, L(t)), with L(t) = g(t)Lo, we change the spatial coordi- 
nate x = yL Q /L(t), where L = L(to), in order to map the problem into the interval (0, L Q ). 
This transformation counterbalances advection, and so the resulting equation reads 



dh ( L \ 2 ^d 2 h Q , „ L _ . 

* = {m) D w-t +F+ ^ x ' t) - (7) 

where we have used the fact that the noise is delta correlated. The dilution term has become 
foV-v = —{g/g)h. It has been disregarded in reaction-diffusion systems due to its irrelevance 
in this context [121 ]. but we will keep it here, where it will show its measurable effects 
on the dynamics. Indeed, dilution has a transparent physical meaning: as the substrate 
grows the deposited material becomes distributed in a larger (d— dimensional) area. This 
matter redistribution causes in turn the propagation of correlations, additionally to diffusion, 
resulting in a different dynamical scenario as the following sections will show. Now we assume 
that the growth function adopts the power law form g(t) = (t/to) 1 , where the growth index 
7 > 0, to find 

The growth index 7 is a new degree of freedom of this problem; it cannot be deduced from 
the other model parameters, and has to be measured directly from the physical system under 
study. Our next step is to assume no flux boundary conditions d x h(0,t) = d x h(L ,t) = 0, 
both due to their physical relevance and because they break translation invariance. It will 
be interesting to see how translation invariance is recovered as a consequence of decorrela- 
tion for rapidly growing domains. Lets decompose the solution in the basis formed by the 
eigenfunctions of the Laplacian on the domain under consideration 

h(x,t) = YX(0cosf^y (9) 

to reduce Eq. ([8]) to a stochastic differential equation for the different modes 
dh n n 2 -n 2 ( t \ 27 „ , 7 , f t x " 



dt ~ L\ W J """" f" n ' V t 



^ Dh n -ih n +[^\ £«(*), (10) 



if n 7^ 0, and 



dho 7 A 



In these equations £ n (t) is a Gaussian random variable with zero mean and correlation given 
by 

(U(t)Uf)) = ^S mn 5(t-t'), if n,m^0, (12) 

(6>(t)£n(0> = 0, if « ^ 0, (13) 

(&(*)&(0> = f *(*-0- (14) 

^0 
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III. EDWARDS- WILKINSON DYNAMICS 



We can straightforwardly derive the equation of motion for (ho) 



dt t 

whose long time solution reads 



~--T{ h o) + F, (15) 



(h ) = ^—t. (16) 
7 + 1 



For the second moment we find 



and the corresponding long time solution 

where the second summand in the right hand side of this equation will be explicitly sup- 
pressed due to its subleading character, but it will be taken into account implicitly later on 
in order to construct Dirac delta functions out of infinite series. For the other modes we 
find 



d(h n ) > > 2 - > 

dt 

and we integrate it to find, in the long time limit, (h n ) — > 0. The correlation obeys the 
equation 



^ ( h m h n ) 



' (m 2 + nyff 2 1 
Llt*< t 



2et< 

(h m h n ) + -r^-Smn, (20) 

Lot ' 
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which solution is readily computable in terms of Misra functions [13|; however its concrete 
form is complicated and not particularly illuminating, so we will not reproduce it here. The 
long time asymptotics of this solution depends on the value of 7. For 7 < 1/2 the interface 
dynamics is dominated by diffusion, and one finds 

(w= Ji m : L i A f) 1 ■ (21) 

D[m l + n 2 j7i z \t Q J 

If 7 > 1/2, then the interface dynamics is dominated by dilution, and the asymptotic solution 
reads 

<^> = (frfe(£f- < 22 > 

Now we can reconstruct the first moments of the solution in coordinate space. For the long 
time mean value we have 

(h(x,t)) = ^—t, (23) 
7 + 1 

and for the long time correlation when x 7^ x' 

= g_ cos (__) cos (__), (24) 

if 7 < 1/2, and 

if 7 > 1/2 and asymptotically for long times, where we have used the decomposition of the 
Dirac delta function in the same basis as in Eq. and we have implicitly taken into account 
the subleading term in Eq. (|18p . Changing back to the original Lagrangian coordinates, 
y = xL(t)/L , we find 

<%, t)h(y>, t)) = + -^-My - y'\ (26) 

(7 + \y 7 + 1 

the solution reduces to random deposition for long times. We thus see that the surface stays 
uncorrelated if 7 > 1/2, while it becomes correlated when 7 < 1/2. We can further analyze 
the correlated phase summing Eq. (1241) to find 

F 2 .„ e 



(h(x,t)h(x',t)) = . t 2 + 



(7 + 1) 2 12DL 



-(—*) [2L 2 -QL max(x,x) +3(x 2 + x' 2 )] , (27) 
\*o/ 



and in y coordinates 

(h(y, t)h(y', t)) = J^f t2 + l2 D L {t) ^ 2L{t)2 ~ 6L(t) max(z/ ' V>) + ^ + ' (28) 
The one point correlation function reads 

<*(*■ m = (^ip' 2 + j^r (Q 1 w - + fe 2 ] • (29) 



or alternatively 



tf) = + [L{tf - 3L(t)y + 3y 2 } . (30) 



Finally, we can calculate the height difference correlation 

([h(y,t)-h(y',t)] 2 ) = ^ D \y-y>\, (31) 

in agreement with the long time classical EW equation. Note that the one point correlation 
function is not spatially homogeneous due to the no-flux boundary conditions (periodic 
boundary conditions or unbounded domains preserve the spatial homogeneity); to find a 
coordinate independent result we define the surface width W(t) as in 

1 f L ° r „ , ,o , e ft V . 



W(t) 2 = -j o [(Kx 1 t) 3 )-(h(x,t)) v \dx = ^ 5 \^-j Lo = ^L(t). (32) 

From these formulas one finds that the width grows without achieving saturation with an 
exponent floo = 7/2, and its square depends linearly on the substrate size, what allows 
to define the roughness exponent a = 1/2, as in the non-growing substrate case, if we 
accept this definition despite the absence of saturation. The marginal situation 7 = 1/2 is 
characterized by an equilibrium of diffusion and dilution, as can be seen by regarding Eq. 



( 120]) . In this case one can see that the correlation length in the x coordinates is A ~ ^/Dto, 
what yields a correlation length in the y coordinates A « \f~Dt. This shows that the fraction 
of interface that becomes correlated is 



A(t) _ VDto 

revealing that the interface becomes globally correlated only if diffusion is large enough, 
or alternatively if the initial system size and growth rate are small enough; otherwise, the 
interface only becomes partially correlated. 
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It is clear that Eq. (1221) is valid if the system is observed from spatial distances \x — 
x'\ ^> (t/to) 1 ^ 2 ^ 7 - This might constitute a good approximation for the two point correlation 
function, but it definitely breaks down when considering the one point correlation. In 
this case we have to consider again the solution of Eq. (TSUI) , but this time in the range 
n,m > (t/t ) 7-1 / 2 - In this scale the interface is again dominated by diffusion and dilution 
may be disregarded, and we obtain the same long time solution as in Eq. (|2T|) . The 
asymptotic behavior is obtained summing this expression, but including in the sum only the 
modes that contribute to the short scale behavior n,m > (t/t ) 7-1 / 2 . The result is 

in the long time limit, where we have used the asymptotic expansion of the Euler-Maclaurin 
formula representation of the series. This result is valid as long as 7^ x 7^ L , otherwise 

(MO, t) 2 ) - (h(0, t)) 2 = (h(L , t) 2 ) - (h(L , t)) 2 « ^ (±) V2 . (35) 

Note that the behavior is now perfectly homogeneous, due to the uncorrelated character of 
the interface in this case. The only exception are the boundary points x = 0, L , because 
they are affected by growth only along one direction, and as a consequence their auto- 
correlation is twice the value of the auto-correlation of any other point not located at the 
boundary. These results hold independently of the reference frame, be it Lagrangian or 
Eulerian, precisely because the interface is uncorrelated. This allows us calculating the 
height difference correlation function 

M-WOf)^^, (36) 

when y 7^ y', 7^ y 7^ L(t) and 7^ y' 7^ L(t). This shows the agreement with the short time 
classical EW equation. Note that the dependence on the system size (both initial and time 
dependent) is the same as in the static domain case, allowing the definition of the roughness 
exponent, which takes on its classical value. 

In order to clarify the 7 > 1/2 dynamics more we will calculate the two point correlation 
function in the long time limit and short spatial scale \x — x'\ (t/to) 1 ^ 2 ~ 7 : 

7 



(h(x, t)h(x, t)} — (h(x, t)) 2 = \x — x'\ ( — ) Tew 



\x — X 



t N ^ 1/2 



(37) 
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or alternatively 



(h(y,t)h(y',t))-(h(y,t)) 2 
where the scaling function reads 



D Z7T Z U 



y'\F< 



EW 



\y-y 



to 



-1/2 



(38) 



U7r\ 

To) 



2tt 



-si 



U7r\ 

L~o) 



(39) 



where si(x) = — j°°[sxn(s)/s]ds is the sine integral. These results are again reminiscent of 
the scaling behavior of the classical EW equation. It is remarkable how translation invariance 
appeared in these formulas, despite the presence of boundary conditions, as a consequence 
of decorrelation. 



IV. GENERAL LINEAR LANGEVIN EQUATION 

We now move to a more general situation in which we consider an arbitrary diffusion 
operator of order ( and an arbitrary spatial dimension d, see Eq. (jSJ). From now on the 
d— dimensional coordinates will be denoted x — > x and y — > y for simplicity. In this case, 
we can proceed exactly in the same way as in the one-dimensional situation to find, instead 
of Eq. (jHJ), the equation 

d t h = d (^j ' 7 \v\ c h -Y h + F+ (j) ^ ^ *)> ( 4 °) 

where the fractional operator |V|^ acts in Fourier space as 

(I = -^An, (41) 

L o 

where (•)„ denotes the corresponding Fourier transformed quantity, and n = (m, • • • , no). 
However, the solution of a boundary value problem for an arbitrary fractional power of 
the Laplacian is obtained not so straightforwardly. For some values of the parameter ( 
(such as for instance ( £ (1,2)) this operator describes Levy flights dispersal, for which 
even the definition of boundary presents difficulties [14]. To overcome this pitfall we will 
decompose the solution using the same basis as in the EW equation case, see Eq. (jSj), in 
the case of a d— dimensional cubic box. This harmonic decomposition implies the no flux 
boundary conditions not only for ( = 2, but for all positive even integer values of this 
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exponent. Indeed, for these values of £ it is easy to check that expanding the solution in the 
d— dimensional analog of (Q is equivalent to prescribe the vanishing of all odd order smaller 
than ( derivatives of the solution at the boundary. This fact allows us to propose a plausible 
definition of solution to the no flux initial-boundary value problem (HOT) as the solution to this 
differential equation which is expressed in terms of the d— dimensional version of expression 
(Q. This way the solution to the initial-boundary value problem is defined in the context 
of evolution semigroup theory as an initial value problem in a prescribed functional domain 
fis| . which in this case is generated by the selected eigenfunctions of the Laplacian, as in 
Eq. ©. 

In the present case the superuniversal threshold above which the interface becomes un- 
corrected turns out to be 7 = For 7 greater than this value (and for large spatial scales 
in case of the second moment), the two first moments of the function height are given by 

<%.«» = s^T*. < 42 > 
W*«W-*» = (^V + S^T*-^ (43) 

Lets now take a look of the one point correlation function. We can proceed in the same way 
as in the previous section to find 

(h( X ,t)*) - (h( X ,t))* ~ L<- d m £ u<~ L °~ d f ' (44) 

by means of the asymptotic expansion of Euler-Maclaurin formula for the series in the last 
step, and where 

E=E-E- («) 

n ni 

In Eq. (1441) we have assumed x = for ease of analysis, as we know that, in the limit 
of uncorrelated interface, all points are equivalent up to a numerical prefactor present in 
the boundary points the rough interface inequality ( > d was assumed too. When 
this inequality is reversed ( < d, then the series does not converge. The situation is the 
same in the case of a nongrowing domain, where for £ < d the interface is flat or at most 
logarithmically rough, but this series is divergent. Note that the one point correlation Eq. 
OH]) grows as t 213 , where j3 is the corresponding growth exponent of the model without 
domain growth, and as Lq q (and also as L(t) 2a ) for the corresponding classical roughness 
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exponent a. These exponents arise as a consequence of the local behavior of the model, 
which is analogous to the static unbounded domain one, see below. 

The behavior of the one point correlation function also helps establishing some dynamical 
properties of the interface when 7 < l/( 

/ f \ (C-dh 

(h(y,t) 2 ) - (h(y,t)) 2 ~ {-) . (46) 

As in previous cases, we note that there is no saturation, unlike in the non-growing domain 
situation, for any 7 > 0. One can define the roughness exponent a from the second moment 
dependence on the system size for long times. In the case of a growing domain one would 
have in principle two possible choices: the initial system size L and the time dependent size 
L(t). It turns out that both yield the same value a = (£ — d)/2, what allows an unambiguous 
definition of this exponent. It is worthy noting that this exponent is the exactly the same 
as in the regime with 7 = 0, with the difference that this latter case is characterized by the 
saturation of the fluctuations. If we allow the definition of the long time growth exponent 
Poo as the power law dependence of the second moment on the temporal variable for long 
times, we find (5^ = r y(( — d)/2. To calculate an effective dynamic exponent z e s we note that 
the correlation length A travels as A(t) ~ (Dt) 1 ^, and so the correlated interface fraction 



at time t is 



A(t) ^ D^t 



L{t) L 

what leaves us 



t 1 ^, (47) 



^eff = ~j ~T ? (48) 

if 7 < 1/C and 00 if 7 > One sees z e &{^ = 0) = ( and we recover the classical case, and 
lim^i/^- z c s(j) = 00, showing the limit in which the interface becomes uncorrelated. 

In the marginal situation characterized by 7 = l/( diffusion and dilution balance each 
other and so the resulting dynamics is given by the concrete values of the equation parame- 
ters. As we have seen, the effective dynamic exponent becomes divergent, and the resulting 
fraction of correlated interface is 

$}"^ 

what shows that for large diffusion and small initial system size and growth rate the interface 
becomes globally correlated. Alternatively, for small diffusion and large initial system size 
and growth rate the interface becomes only partially correlated. 
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As we have seen, the roughness exponent can be defined a = (£ — d)/2 independently 
of the value of 7, as in the static domain case. For a supercritical 7, the long time growth 
exponent reads /3oo = (1 — d/Q/2, equalling the growth exponent in the classical case. For 
any 7 7^ 1/(2C) the inequality a 7^ PooZes holds, while the equality a = (3^ is true only 
if 7 > On the other hand, we have the alternative relation a = /^oo-^oo, for z^ = 

max{l/7,£}. Indeed, there is a close connection among I/7 and the dynamic exponent, i. 
e., this quantity describes quantitatively the speed at which correlations propagate along 
the interface. This relation is further explored in the following sections. 

One can clarify things further by calculating the two point correlation function for those 
bulk points that lie closer than the correlation length \x — x'\ <C (Dt) 1 ^~' y . In d = 1 we 
obtain 



(h(x,t)h(x',t)) - (h(x,t)} 2 
or alternatively 

(h(y,t)h(y',t))-(h(y,t)) 2 = \,j ,j 
where the scaling function reads 







c-i 






\x — x'\ 


(a; 


T 


\x — x'\ 


(a 













7-1/c 



(50) 



/|C-i 



\y -y \ \ — 



-i/c 



(51) 



eL 



C-l 

u l -< 



4tt<D 



cos 



U1TS 



ds, 



(52) 



which integral can be considered as a trigonometric variant of the Misra function [l3]. These 
formulas have been found in the long time limit after adiabatic elimination of highly oscilla- 
tory functions (which results from a direct application of the Riemann-Lebesgue lemma). Re- 
sult (jHJ) together with (15T|) allows us to recover the classical critical exponents a = (( — l)/2, 
(3 = 1/2 — l/(2£), and z = (, when we consider Lagrangian coordinates y and distances 
shorter than the correlation length. Note that Eqs. (150]) and (|5T|) only depend on \x — x'\ 
and \y — y'\ respectively and thus they are translational invariant: this is a consequence of 
decorrelation, what makes the interface bulk behave similarly to the case of an unbounded 
domain for short spatial scales. In this limit, boundary conditions are not affecting the 
dynamics of bulk points, but boundary points show a different prefactor as in the strictly 
local situation, see (135!) and 
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V. TEMPORAL CORRELATIONS 



In order to calculate the temporal correlations we need to consider EW dynamics Eq.® 
in the short time limit, where the growth exponent j3 becomes apparent. The homogeneous 
solution of its Fourier transformed representation Eq. lflU]) is 



t 



h n(t) = ( — ) exp 



-7 



n 



2 TT 2 D t^t 1 '^ - tr 



h n (t ) = G n {t)h n {t ), (53) 



L 2 I-27 

that yields the following complete solution when the initial condition vanishes: 

•t r+ \ 7/2 



Kit) = G n (t) J G~\t) (^j 1 Ur)dr. (54) 
The one point two times correlation function then reads 

r, /•min(t/) / . \ 7 

(h n (t)h n (t')) = J-G n (t)G n (t') jf G-\t) l^j dr, (55) 

and after inverting Fourier we arrive at the real space expression 

(h(x,t)h(x,t')) = J2(hn(t)K(t')) cos 2 . (56) 

The propagator G n (t) suggests the scaling variable v n ~ nt 1 / 2-7 in Fourier space, that 
corresponds to the real space scaling variable u ~ xt y ~ 1 / 2 , as can be read directly from 
the last equation. This again suggests the definition of the effective dynamical exponent 
z e ff = 2/(1 — 27). If we express the correlation Eq. ( J55l) for t — t' in terms of the scaling 
variable v n (and we refer to it as C{v n ) multiplied by a suitable power of t) and we introduce 
the "differential" 1 = An ~ t 7-1 / 2 At>, we can cast the last expression in the integral form 

(h(x, t) 2 ) - (h(x, t)) 2 = t 1 ' 2 ^ C(v n ) cos 2 (j^j dv n , (57) 

where the series converges as a Riemann sum to the above integral when 

7- 

l 



Dt<(Ll + Dto)-^, (58) 

tn 



or equivalently t -C t c ~ Lq" s , for t c being the time it takes the correlations reaching the 
substrate boundaries, assuming that the substrate initial size is very large. If 7 < 1/2, 
the whole substrate becomes correlated, yielding a finite t c ; for 7 > 1/2 the convergence 
of the Riemann sum to the integral is assured for all times, corresponding to the physical 
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fact that the substrate never becomes correlated. In front of the integral we find the factor 
t 1//2 , compatible with the growth exponent (3 = 1/4, and the integral can be shown to be 
absolutely convergent due to the Gaussian dependence of G n (t) on n. 

The general situation in which we deal with a d— dimensional substrate and the diffusion 
is mediated by an operator of order ( can be constructed along the same steps. In this case 
the propagator reads 



G n (t) = [ 1 



— dq 



exp 



n 



<<K<D t^t 1 ^ - in 



7C 



(59) 



suggesting that the scaling variables are v n ~ ni 1 '^ -7 and u ~ xt 7 ^ 1 ^, and a definition 
for the effective dynamical exponent z e ^ = £/(l — 7^). We find again convergence of the 
Riemann sum to an integral for any time if 7 > and for short times 

Dt < (Lo + Dio)^, 

if 7 < 1/C, in agreement with the expression for the correlation time t c ~ L Z q s . This integral 
is again absolutely convergent as it decays superexponentially for large values of the scaling 
variable v n , leading to the result 



(h(x,t) 2 ) - {h{x,t)f 



t 



i-d/c 



(60) 



in agreement with the classical growth exponent (3 = 1/2 — d/(2(). 
We are now in position to calculate the temporal auto-correlation 



A(t,t') = 



{h(x,t)h(x,t')) 



min{i, t'} 



(61) 



(h(x, t) 2 )l /2 (h(x, t' ffj 2 Vmax{i, t'} , 

where A is the auto-correlation exponent and (-)n denotes the average with the zeroth 
mode contribution suppressed, as in ( 1601) . The remaining ingredient is the correlation 
(h(x,t)h(x,t'))o. Going back to Eq.( |56l) we see that the Fourier space scaling variable now 
reads 

Vn= V 7 U - 62 

1 - 7( 

If 7 < l/( the term max{i, i'} 1-7 ^ is dominant and the factor in front of the convergent 
Riemann sum reads 

max{t,t'}~ d/<: imii{t,t'}, (63) 
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after the time integration has been performed and in the limit max{t,t'} ^> min{t, t'}. In 
this same limit, but when 7 > the term min{t, t'} 1 " 7 ^ becomes dominant and the 
pre-factor reads 

max{t, t'Y dl min{t, ^}i-*/C+*r. (64) 

The resulting temporal correlation adopts the form indicated in the right hand side of (I6TI) . 
where 

xJ" + dK " 7<1/C («) 
\fl + ld if 7 >1/C, 

or alternatively 

A = /?+-, (66) 
where (3=1/2 — d/(2Q and the new dynamical exponent is defined as 

z x = min{C, 1/7}- (67) 



This last form is the natural generalization of the corresponding one in 17|, and tells us 
that correlations are propagated either by diffusion or dilution: the dominant mechanism is 
chosen in each regime. We can extract more information about the correlation function, as 
it decay properties 

(h(x,t)h(x,t')} ~ max{t,t'}- d/z , when max{t, t'} -> 00, (68) 

signaling that it decays to zero for long times. Also, the short time behavior of the auto- 
correlation function is 

/ mini* t'\ \ 1 ~ d/C 
A{t, 0-1-^1- \ { , when max{t, t'} w min{t, if}, (69) 



max{t, t'} 

homogeneously in 7, where R = R(£,d,j) is a universal function of its arguments. This 
behavior is compatible to the one found in the 7 = case [3]. It indicates that the short 
time properties of the auto-correlation are independent of the substrate growth velocity, but 
the long time behavior is influenced by the mechanism by which correlations are propagated, 
be it diffusion or dilution. 



VI. PERSISTENCE 



The persistence of a stochastic process denotes its tendency to continue in its current 
state. When considering the dynamics of a fluctuating interface, one refers to the persistence 
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probability P + (ti,t 2 ) (P-(ti,t 2 )) as the pointwise probability that the interface remains 



13, 



181 ] . Herein, as in [8J], we concentrate 



above (below) its profile at t% up to time £2 > £i 
on the case in which the initial profile is flat, and we suppress the contribution coming 
from the zeroth mode as in the last section. For the stochastic growth equations under 
consideration the symmetry h n — > — h n for all Fourier modes n > holds, implying the 



equality P + = P_ = P. For long times £2 3> t\ we have the power law behavior 17|, [18] 



P(*i,t 2 ) ~ (ti/t 2 ) e , (70) 

defining the persistence exponent 9. It was previously calculated in the limit ( — > 00 when 
7 = |is| 

up to higher order terms, and in this same limit when d — 1 and 7 = 1 

« - - (72) 
2 2( 

up to higher order terms. The goal of this section is to calculate the persistence exponent 9 
in the limit ( — > 00 for a finite, but otherwise arbitrary, value of 7. In order to proceed with 
the calculation, we need to consider again the normalized auto-correlation function, i. e., 
the left hand side of (I6"TT) . This time we will not focus on the limit max{t,t'} ^> mirvfi, t'}, 
instead we will consider an arbitrary relation among t and t'. In this case we have [l9( 



(h(x,t)h(x,t')) ~ max{t,t'}- d7 min{t,t'} (max{t, t'} 1 ^ 2 - min{t, t'} 1 ^) 



/\l--yz\ —d/z 



min{t, t'} maxjt, t'} 



d/z 



2 



ryd—ld jd—1 min{£, t 1 } max{t, t'} 



JZ 



(73) 



.7-2-1 z jz-r t{t')^ z + v z t' 

where 2 Pi(^i, x 2 ; x 3 ; x 4 ) is Gauss hypergeometric function [2^]. In order to derive the per- 
sistence exponent we consider the auto-correlation function in logarithmic time T = ln(t) 

A(T,T>) = ^^tln^^", ~ e-^ T ~ T \ when \T-T'\ -> 00. (74) 
(h(x,eT) 2 )l /2 (h(x,en 2 )l /2 

Note that this is the correlation function for the normalized (to unit variance) function 

height, which becomes stationary in the logarithmic temporal variable. For short times we 

have 

A(T, T') — l + O (|T — T'| 2/3 ) , when |T-T'|->0, (75) 
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homogeneously in 7. The first order term is a power 2(3 = 1 — d/( < 1, classifying 



the process as a Slepian non-smooth one 21(. This fact together with the asymptotics 
A(T,T') ~ e -( 1 / 2 + d i)\ T - T '\ when ( — > 00 means that we can calculate the persistence expo- 
nent 9 perturbatively about 9 = 1/2 + dj for large ( (and finite d and 7) in the following 



fashion 



9 



3 



1 



+^7 



1 - 



1 + 2d 7 



7T 



(A(r) - 



-(l/2+d 7 )r 



} |! _ e -(l+2d 7 )r} 



-3/2 



(76) 



for r = |T — T'|, yielding, in the limit £ — ► 00, the result 





1 , d 

— h 07 . 

2 2C 



(77) 



up to higher order terms. This last result is reminiscent of the one obtained in 8J, but it is 
corrected by the effect of dilution, as we shall discuss in the next section. We can see that, 
in the limit considered, the interface is less persistent than in the case of a static domain. 
This is so because in the uncorrelated phase dilution acts as a relaxation mechanism on the 
strictly local scale with a higher efficiency than diffusion. Note that in the limit 7 — > one 
does not recover the static domain result [3] . This is so because in this calculation we have 
assumed 7 > l/£, and so the vanishing 7 limit implicitly implies a faster vanishing l/( limit. 
As a conclusion we find that when 7 — > the persistence exponent 9 — > 1/2. 

VII. CONNECTION TO RADIAL GROWTH 



As mentioned in the Introduction, one of the characteristics of radial growth is its growing 
domain interface. Herein, we will use the results derived in previous sections to put former 

nn 

derivations in radial geometry [4], [8J in a broader context. In this case as well dilution plays 
an important role on the interface dynamics. As radial interfaces grow in time, the interfacial 
matter becomes diluted among the new deposited matter and the correlations transported 
simultaneously. Note that the physical origin of dilution here is the same as in sec. HIl and 
thus it is not related to the surface curvature. 

The one- dimensional radial counterpart of the general linear Langevin equation (|40j) could 
be defined as jg] 

= jFt^ 1 + -\V e \ C r + -^ V (9,t), (78) 
at r<> Jr 
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for the field r(9,t), where reparametrization invariance [3] has been taken into account, but 
dilution has been disregarded. Its analysis yields, for 7 > l/£ and performing a van Kampen 
system size expansion about the homogeneously growing state, the long time large angular 



scale correlation function 



23] 



for 7 < 1, 
for 7 = 1, and 



(r(6,t)r(8',t)) ~ t 1 ^ 6(6 - 9') ~t5(s- s'), 



(r(9, t)r(9', t)) ~ ln(t) 5(9 - 9') ~ t \n(t) 5(s 



(r(9, t)r(9', t)) ~ t 1 ^ 5(9 - 9') ~ f 5(s - s'), 



(79) 



(80) 



(81) 



for 7 > 1, where s — s' ~ P(9 — #') is the arc-length scale. For 7 < 1 we recover the random 
deposition correlation, while for 7 > 1 we found, in the arc-length variable, an average rapid 
roughening version of it, this is: 

if 7 < 1, 



s(6»=2tt) 



(r(s,t)r(s',t)) ds ~ < 



s(0=O) 



t ln(t) 



(82) 



if 7 = 1, 
if 7 > 1. 

This result emerges when the dilution term is not taken into account. When we contemplate 
the effect of dilution, as in the previous sections, we find a pure random deposition correlation 
in Lagrangian coordinates 

(r(s,t)r(s',t)) ~t5(s-s'), (83) 

homogeneously in 7 (provided 7 > 1/C) ; as in (H3]) . Correspondingly, the prefactor of the 
Dirac delta is t 1-7 in the Eulerian setting. One can see that the large scale results derived 



for radial interfaces 



231 ] are identical to the ones found here for homogeneously growing 



domains, once dilution is introduced. These results can be straightforwardly generalized to 
an arbitrary dimension d. In this case one needs d angles to parameterize the interface in the 
Eulerian setting, that will lead to d different arc-lengths in Lagrangian coordinates and to 
a d— dimensional Dirac delta specifying the spatial properties of the uncorrelated interface. 
As a consequence one finds 



(r(s,t)r(s',t)) ds 



t 

t ln(t) 



if 
if 
if 



7 < 1/d, 
7 = 1/d, 
7 > 1/d, 



(84) 
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when dilution is not considered; s = (si, • • • , so) is the set of all arc-lengths and the integral 
extends to the whole domain. If we contemplate dilution, then the resulting correlation in 
the large spatial scale and for long times is pure d— dimensional random deposition. The 
physical reason for this enhanced stochasticity when dilution is not present is the following. 
Dilution distributes the fluctuations along the interface at the growth rate, keeping the total 
amount of noise constant. In its absence, the interface is composed of unconnected sites, 
and new ones are added in the process of domain growth. They act as independent sources 
of noise, and so they contribute to augment the overall fluctuations. 

The study of radial growth performed in 8j on the strictly local scale does not contemplate 
the effect of dilution either. For the shake of completeness, and to facilitate comparisons, 
we have derived the previous sections results deliberately disregarding the effect of dilution. 
As expected, the long time short spatial scale behavior is reminiscent of the unbounded 
substrate situation, as happened in ([50]) and (15T|) . On the contrary, the large spatial scale 
dynamics shows a different phenomenology characterized by a particular random deposition 
effective behavior. Equivalently, the long time interval (max{i, t'} 3> min{t, £'}) asymptotics 
is affected by substrate growth, and the absence of dilution modifies the corresponding 
results. If 7 < the correlation function decays to zero at infinity as a power law 

(h(x,t)h(x,t')) ~ (max{t,t'}) d7-d/c , when max{t, t'} -> oo, (85) 

while for 7 > l/£ this correlation function approaches a non-zero value. This fact is related 
to the behavior of the Fourier modes (1591) . which in the absence of dilution decay to zero 
for 7 < 1/Q in the long time limit, while for 7 > l/£ they approach a non-zero value 
asymptotically in time, analogously to the particular case analyzed in Q|. 
The temporal auto-correlation function, in the absence of dilution, is 

(h(x, t)h(x, t')) ( min{t, t'} x A 



where 



(h(x, ml /2 (h(x, t') 2 )l /2 Vmax{t, f} 



+ W~dl if 7 < 1/C, 



2 1 2C 

1 d_ 

2 2£ 



' " if 7>1/C, 



(86) 



(87) 



or alternatively 



A = /3 + —, (88) 

Zett 
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where = 1/2- d/{2Q and 

JC/(1-7C) if 7<1/C, 

^eflf = < (89) 

I oo if 7 > 1/C- 

From these formulas one can clearly read that when dilution is suppressed there is no mech- 
anism for correlation propagation and thus the system behaves as an effective particular 
random deposition model in this limit. Similar information can be obtained from the per- 
sistence exponent, obtained this time perturbatively about 8 = 1/2: 

d ~\-^ C - °°> ( 9 °) 



generalizing the previous result [8j. 
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'his, together with numerical simulations suggesting 



9 ~ 1/2 almost homogeneously in ( [8J|, reinforce the idea of an effective particular random 
deposition behavior. However, for ( large enough 9 < 1/2, and in consequence the process 
is more persistent than random deposition. This fact admits a transparent physical expla- 



nation. For a static domain 9 decreases for increasing (3 13]: the reason is that in this case 
the exponent f3 contains information on the relaxation properties of the interface (through 
its dependence on the dynamic exponent z). In this case relaxation is mediated by diffusion, 
that connects the different interface points and thus pushes the interface towards its mean 
value, diminishing the persistence of the fluctuations. For a growing domain, dilution acts as 
the relaxation mechanism when diffusion becomes inoperative (in the uncorrelated phase). 
Suppressing dilution, there is no relaxation mechanism left, following that the exponent (3 
only contains information about the strictly local fluctuational properties of the interface. 
Indeed, the interface variance grows as (see sec. lIVp . and so for smaller j3 we have 
weaker fluctuations intensity, implying a longer first passage time. This implies in turn a 
smaller value for the persistence exponent 9, as persistence is nothing but a first passage 



problem 181 ] . This explains how 9 may increase for increasing /3, although understanding 



the whole numerical sequence of values for the persistent exponent in [8| would require a 
deeper analysis. 

These last results for the temporal auto-correlation and persistence show that in the un- 
correlated phase dilution is the only responsible for correlations propagation. Interestingly, 
the simulations performed in [g] for the Eden model show a temporal autocorrelation 
function fully compatible with the one described here for an uncorrelated interface in the 
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absence of dilution. The two persistence exponents, for above and below the mean fluctua- 
tions (which are different for the less symmetric Eden interface), are greater than 9 = 1/2, 
but are notably smaller than the static domain ones 8J. This result is surprising because for 
a radial cluster grown according to the Eden rules, to which new cells are added at random 
positions on its interface, dilution is expected to occur. It would be very interesting to un- 
veil the mechanism counterbalancing dilution in this case. It might be a consequence of the 
particular way in which an Eden cluster grows, or perhaps due to some possible nonlinear 
effect acting on the interface as a consequence of geometry. 



VIII. CONCLUSIONS 



In this work the dynamics of linear stochastic growth equations whose domain size grows 
in time as power law t 7 has been studied. The growth index possesses one critical value 
7 = 1/C, for ( being the order of the diffusion operator. If 7 < l/( the interface correlations 
are propagated by means of diffusion at a faster speed than domain growth, resulting in a 
fully correlated interface. The time it takes correlations to travel the whole interface depends 
on the initial substrate size t r ~ Lq s for z c g — C/(l — tC)j or alternatively t c ~ L^t^^. 
The roughness exponent a can be defined from the strictly local properties of the interface 
uniformly in the growth index. This value of a is exactly the same one that is obtained 
in the static domain case. For any 7 > saturation never occurs, and the interface width 
continues to grow for all times with a long time growth exponent /3oo = 07 when 7 < 
note that in this regime /3oo < (3 so there is partial saturation of the fluctuations. The relation 
Ax> = «7 also implies that correlations travel like t 1 in the long time regime, after they have 
spread globally on the interface. Prior to that they propagate as t 1 ^, because diffusion is 
a faster mechanism for information transfer. Once they have reached the interface limits, 
this transfer speed is limited by the slower process of domain growth, resulting the stated 
exponents relation. This phenomenology is independent of whether we contemplate dilution 
or not: it is a strictly direct consequence of domain growth, not of dilution (although dilution 
carries on information at the same velocity). 

The regime in which 7 > l/( is characterized by a loss of correlation along the interface. 
This translates into a delta correlated spatial correlation for long times and large spatial 
scales. The correlations for short spatial scales and time intervals are reminiscent of the 
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ones found in the case of a static unbounded domain, revealing that diffusion is acting 
at this level. On the other hand, large spatial scales and long time intervals both reveal 
that the dominant mechanism for correlations propagation is the dilution effect created by 
domain growth. By means of dilution, correlations travel at the same speed at which the 
domain grows, so a global correlation of the interface becomes impossible. The situation is 
further clarified by the calculation of the auto-correlation and persistence exponents. The 
autocorrelation exponent A = (3 + d/( for all 7 < 1/C (including 7 = 0) shows that one 
site interacts with itself at former times by means of the growth process (indicated by the 
first summand (3) and with neighboring sites by means diffusion (indicated by the second 
summand d/(), which is dominant in this regime. If 7 > then the auto-correlation 
exponent reads A = (3 + d*y. This illustrates how, for fast domain growth, dilution replaces 
diffusion and becomes responsible for the interaction with the neighboring sites. A similar 
conclusion is reached by analyzing the persistence of the surface fluctuations. While the 
auto-correlation exponent yields information about the long time interval asymptotics, the 
persistence exponent carries on complimentary information obtained from averaging over all 
possible time interval lengths. It reads 9 = l/2 + d r y — d/(2() when ( — > 00 and 7 > 1/C, and 
so it is greater than in the static domain situation, what implies that the interface is less 
persistent. A reduced persistence is associated with a stronger tendency to go back to the 
mean, which is mediated by a stronger coupling with the neighboring points through dilution, 
more efficient than diffusion in the limit considered. Additionally, the persistence exponent 
increases with the growth exponent: a higher (3 corresponds to stronger fluctuations and 
thus to a shorter first passage time. This is in contrast to what happens in the 7 = 
situation, where the persistence exponent decreases for increasing (3 [17|. The reason for 
this is that in this case f3 contains information about the strength of the coupling among the 
interface sites. Smaller (3 implies a stronger coupling and a correspondingly less persistent 
interface. 

The crossover situation 7 = 1/C is characterized by a nonuniversal behavior. There is a 
strong dependence on the parameter values that enter in competition to yield the resulting 
system dynamics. In the generic case 77^ 1/C, a number of the results derived in this paper 
can be connected to the Family-Vicsek ansatz Q by means of the simple substitution 
L — > L(t) taking into account the temporal evolution of the system size directly in this 
ansatz. From the results of this paper, one would expect this to be so in the regime 7 < 1/C, 
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or for 7 > 1/C provided short spatiotemporal distances are under consideration. When 
7 > 1/C, and for long time intervals and/or large spatial scales one would perhaps expect 
a different result. Indeed, in this regime correlations are propagated by means of dilution, 
which becomes more effective than diffusion. However, the surface dynamics is still well 
described by the Family- Vicsek ansatz. To see this, consider for instance the unbounded 
domain situation uniformly growing with a growth index 7. In this case the Family- Vicsek 
ansatz tells us that the two points correlation is: C(\x — x'\,t) = t 2/3 F(\x — s'l/t 1 ^ -7 ) in 
Eulerian coordinates x and 2(3 = 1—d/ ( for the linear stochastic growth equations considered 
here. For long times and 7 > 1/C we find C(\x - x'\,t) = t l -^t d ^- d/c F{\x - x'|/t 1/c " 7 ) 
t l ~ d,y 5(x — x') when t — > 00, i. e., we recover the standard random deposition correlation. 
This result fully agrees with the correlations derived from the stochastic growth equations 
which take into account dilution. On the other hand, the correlations shown at the beginning 
of sec. IVHI were obtained suppressing the dilution term and cannot be derived from the 
Family- Vicsek ansatz. So we see that this ansatz implicitly takes into account dilution. 
This miracle occurs because the Family- Vicsek ansatz neglects the memory with respect 
to the initial condition at to- However, this memory effect is present in the stochastic 
growth equations with no dilution, leading to a different result |22|. Introducing dilution 
asymptotically erases the memory with respect to the initial condition, what implies in turn 
the coincidence of the results from the stochastic growth equations and from the Family- 
Vicsek ansatz. Of course, in this reasoning we have assumed the rough interface inequality 
( > d; otherwise the appearance of non-universal anomalous dimensions is indeed possible 
as we have shown for ( = d in 23]. 

One of the characteristics of radial growth is the growing domain size. We have iso- 
lated this effect, what will hopefully allow a better understanding of the dynamics of radial 
interfaces. In order to facilitate the comparison with previous work we have deliberately 
neglected the dilution term on the dynamics, although this term arises naturally in the sce- 
nario considered. Our results compared favorably and extended those of [sj]. The resulting 
analysis in the fast growth regime showed an interface without relaxation mechanisms, and 
a corresponding temporal auto-correlation function decaying to a non-zero value in the long 
time limit. Once dilution is suppressed, and because diffusion is inoperative in the large 
scale for 7 > there is no remaining coupling among the surface sites, as revealed by the 
auto-correlation exponent A = f3 signalling that one site only interacts with itself at differ- 
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ent temporal points. The information obtained from the persistence exponent is the same. 
Note that the no dilution assumption allows growth characteristics beyond what one would 
expect from the static domain results. For instance, an average rapid roughening effect (not 
real rapid roughening which would imply pointwise estimates out of reach in this context) 
is possible for 7 > 1/d, see Eq. ( I8~4"j) . because fluctuations are no longer being distributed 
along the growing domain. Values of the persistence exponent 9 < 1/2 strictly smaller than 
the random deposition one are possible for large values of (. Both become impossible if the 
domain is non-growing or if we contemplate dilution. Interestingly, simulations performed 
with the Eden model (for which 7 = 1) show that it behaves as if dilution were not present 
8j. Basically, this means that an Eden cluster grows as a random deposition process, with 
all its interface points being completely decorrelated (even at the local scale for long time 
intervals), but with weaker surface fluctuations as given by a smaller value of (3 (— 1/3 for 
this particular model). The only dissenting factor is the numerically measured value of the 
persistence exponents, greater than the expected 9 « 1/2 [8]. Anyway, these values are 
consistent with the rest of results, as they are considerably smaller than the static domain 
ones [nj], implying a more persistent interface. They may be the consequence of a KPZ 
nonlinearity (expected from the [3 = 1/3 exponent) acting on the interface, as it propa- 
gates correlations linearly in time Q| , and so it can compete with the linearly growing Eden 
interface. 

When studying a radial interface, it seems necessary to take into account both 
reparametrization invariance [?| and dilution 12 |. The Eden model seems to be, in princi- 
ple, not an exception, as its interface grows by the addition of new cells randomly placed 
on the cluster surface. However, numerical results pointed to the fact that dilution is not 
operative in this case [sj]. It would be very interesting to understand what mechanism is 
counterbalancing dilution in this model. Some candidates are the presumed KPZ nonlinear- 
ity present in the Eden surface dynamics, which is able to propagate correlations linearly 
in time (exactly the same velocity at which dilution would operate at the linearly growing 
Eden interface), or the nonlinearities implied by reparametrization invariance. 
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